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Abstract. 

We observe a complex change in the hopping exponent value from 1/2 to 1/3 as a 
function of disorder strength and electron density in a sodium-doped silicon MOSFET. 
The disorder was varied by applying a gate voltage and thermally drifting the ions to 
different positions in the oxide. The same gate was then used at low temperature to 
modify the carrier concentration. Magnetoconductivity measurements are compatible 
with a change in transport mechanisms when either the disorder or the electron density 
is modified suggesting a possible transition from a Mott insulator to an Anderson 
insulator in these systems. 
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1. Introduction 

The crossover from Mott to Efros-Shklovskii (ES) variable-range liopping (VRH) lias 
been extensively investigated either to confirm the existence of a Coulomb gap in 2D 
[l] or to experimentally verify its theoretical formalism ^J. It has only been observed in 
the temperature dependence of conductivity of insulating devices close to the metal-to- 
insulator transition. Effectively, in that region, the crossover temperature, a density of 
state dependent entity is a slowly varying function of the electron density and disorder. 
This prevents a change in the value of the hopping exponent as a function of the two 
previous parameters from being observed. However, the existence of impurity band 
tails in disordered insulators [3l HI [5] allows such a study to be performed at a fixed 
temperature and in the insulating regime. In this paper we show that these experimental 
conditions are achieved by using sodium ions in the oxide of a silicon MOSFET. While 
the gate controls the electron density at low temperature, the in-situ modification of 
the position of the ions in the oxide by the simple application of a gate voltage at room 
temperature allows the disorder strength to be changed. We obtained similar results 
to Pepper on small silicon MOSFETs [6] when the ions are close to the interface and 
a comparable electronic situation to Ghosh's GaAs/AlGaAs devices with a fixed silicon 
5-doped layer |7] when they are close to the metal gate. The study is extended to 
intermediate positions by analysing the variation of the source-drain conductivity in 
temperature and in magnetic field. 

2. Experiments 

All measurements were performed on identically processed silicon MOSFETs (Fig. la). 
Such devices have been widely used because of the ability to continuously vary the 
electron density and thus the Fermi energy by use of a metal gate. The geometry of 
the devices was circular (Corbino) to avoid leakage current paths around the source and 
drain contacts. The devices were fabricated using a high resistivity (lO^fi.cm) (100) p- 
silicon wafer to minimize the scattering with boron acceptor impurities, especially close 
to the silicon-oxide interface. A 40- nm thick gate thermal oxide was grown at 950 °C in a 
dry, chlorine-free oxygen atmosphere. Contacts were realized by implanting phosphorous 
at high dose and sputtering aluminium. Sodium ions were introduced onto the oxide 
surface by immersing the device in a 10"^ N solution (- 7 x 10^^ cm ) of high-purity 
sodium chloride in deionized water for 30 s. The surface of the chip was then dried with 
nitrogen gas and an aluminium gate subsequently evaporated. The effective gate length 
of the primary Corbino MOSFET was 2 /xm and the diameter of the interior contact 
was 110 /im. The Na"*" ions were drifted through the Si02 by applying +4 V dc at 65°C 
for lOmins but did not diffuse into silicon|8]. The device was then slowly cooled down 
to about 300 mK. Once the temperature lower than about 150 K, the ions remain frozen 
in position allowing consistency in temperature dependence measurements. A series of 
drifts (Di for the i drift) were performed and allowed modifying the position of the ions 
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in the oxide. The last drift (D5) was carried out at 80°C for Ih. The undrifted device 
(DO) were used as a reference and for estimating the ion concentration in the oxide. 
Standard low- noise lock-in techniques with an amplifier gain of 10^ V/ A were used to 
measure the source to drain conductivity. The ac excitation amplitude was Vac=10/xV 
at a frequency of 23 Hz. The dc offset of the amplifier was suppressed using a blocking 
capacitor. The gate voltage was controlled by a high resolution digital to analogue 
converter and the temperature measured by a calibrated germanium thermometer. All 
experiments were performed in an ^He cryostat where the magnetic field was applied 
perpendicular to the Si-Si02 interface. 

3. Results and discussion 

When the ions are not drifted (e.g. close to the metal gate), the variation of the source- 
drain conductivity with gate voltage is similar to sodium free MOSFETs and, in the 
region of study (subthreshold region), the conductivity depends exponentially on gate 
voltage. Following a drift of the ions, the onset voltage for conduction is shifted towards 
negative gate voltages indicating the presence of positive charges in the oxide (Fig. lb). 
Fluctuations in the conductivity also appear and strengthen as the ions approach the 
Si-Si02 interface. These fluctuations are reproducible both in position and height with 
time and with thermal cycling up to 120 K. This result is in agreement with Pepper's 
results on short channel length silicon MOSFETs with low-doped substrate [6]. The 
absence of an impurity band as seen by many authors [3l |5] may indicate the presence 
of strong potential fluctuations at the interface and a large oxide charge concentration 
often associated with a decrease in the impurity band visibility [9J . Once at the vicinity 
of the interface, any further drift becomes inefficient. By neglecting the ion distribution 
in the oxide, we can estimate the position of the ions for intermediate drifts from the 
experimental threshold voltage shift AV^^ = Vj^ ■ — Vj^Q = e (rf — /ij) A^ox/ (^o e) ? where Aqx 
is the effective oxide charge density, d the oxide thickness and hi and Vj^j respectively 
the distance of the ions to the Si-Si02 interface and the threshold voltage after drift 
Di (Fig. Ic). Threshold voltages were determined by the linear extrapolation to zero 
of a iygj at 290 mK. Allowing for this approximation, we find an upper bound for 
Aqx ~ 1-2 X 10^^ cm~^ for the ion density by supposing that all the ions reached the 
interface after drift D5 (/15 ~0). This gives a corresponding ion separation of d^^ ~9 nm. 
The value for Aqx accounts for the active charges (e.g. not neutralized at the interface) 
and represent about 45 % of the total oxide charge [TO] . It is an overestimate of the 
real ion concentration because it is plausible that some ions have already reached the 
Si-Si02 boundary before performing drift D5. Subsequent drifts would then change the 
distribution of ions at the interface by pulling the remaining ions in the oxide to the 
interface. 



Variation of the hopping exponent in disordered silicon MOSFETs 



4 



Na* 









///// 












Drain 




Source 







0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
-0.1 
-0.2 



p-Si 



Si 








SiO^ 


Al 




-10 10 20 30 40 50 

h (nm) 



-0.6 -0.4 -0.2 0.2 0.4 0.6 



Figure 1. a) Cross-section view of a Corbino MOSFET used in the experiment after 
drift (D5). b) cr (Vg) for sodium ions close to the metal gate (DO) to ions close to 
the Si-Si02 interface (D5) at 290 niK. c) Estimate position of the ions in the oxide for 
different drifts from AVj^, . 




Figure 2. a (T) when the ions were close to the Al-Si02 interface (DO) for Vg = 0.7 V, 
0.57V and 0.5V (from top to bottom), at intermediate position in Si02 after drift D2 
for Vg = 0.6 V, 0.5 V, 0.45 V and 0.35 V (from top to bottom) and close to the Si-SiOz 
interface after drift D4 for Vg = 0.5 V, 0.4 V, 0.3 V and 0.2 V (from top to bottom). 
Optimum hopping exponents are given for information. Best fits are shown with dotted 
lines. 



3.1. Temperature dependence 

We studied the influence of tlie position of tlie ions in tlie oxide on tlie electronic 
properties of the device by measuring the temperature dependence of the conductivity 
after the different drifts. In all cases, the device showed an insulating behaviour in the 
region of study and the conductivity was well fitted to the generalized hopping formula 
over several orders of magnitude between IK and 20 K typically (Fig. 2) [TT] : 

a(T)=aoT-^'^e-(^''/^)^ (1) 
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Figure 3. Variation of s (Vg) after after drifts Di. Dotted lines represent the value of 
the hopping exponent in the Mott regime (s = 1/3) and in the ES regime (s — 1/2). 
The values of s are determined within an error ±0.002. 



Below 1 K the conductivity is mostly determined by scattering with acoustic 

phonons. The fitting procedure used to find the exponent p and s is a standard one based 

on the minimization of the chi-square from the non-Arrhenius plot of ln(crT*'^) versus 

Best fits are obtained with p = 2 in all cases. Such a value for p has already been 

observed in similar devices [3j as well as in Si: As [12j and theoretically predicted [13]. 

When the ions are close to the metal gate, s ~ 1/2, indicating the presence of ES VRH 

[1] (Fig. 3). For intermediate oxide positions, we find a smooth change from s ~ 1/2 (ES 

VRH) to s ~ 1/3 (Mott VRH) when decreasing the gate voltage [H]. After the drifts D2 

and D3, this crossover appears at Vg ~ 0.47 V and 0.45 V respectively. This indicates the 

Coulomb gap progressively diminishes as the electron density is decreased. Effectively, 

the value of T1/2 experimentally increases when the electron density is decreased, such 

/ \ 1/2 

that, at some point, the hopping energy in the ES regime ~ (Tri/2 ) becomes greater 
than the Coulomb gap ~ T^/^/Ti/^ The ES VRH then changes to Mott VRH 

for T > T* ^ (Ep) . Because of the presence of band tails in our device, both 

the density of states at the Fermi energy (-^f) ^^'^ localization length ^ are 
decreasing functions of the electron density, indicating that the ES regime is more likely 
to be present at high electron density. This result is inconsistent with the screening of 
the Coulomb interaction by metallic gate as previously reported in similar devices [31 [17] 
and in other 2D systems where the Mott regime is expected at high electron density |18j . 
After the drift D4, the crossover is still present although 0.32 < s < 0.44. Finally, when 
the ions are close to the Si-Si02 interface (drift D5), hopping transport is still active but 
s (Vgj becomes non-monotonic (0.28 < s < 0.32) due to the presence of fiuctuations in 
the conductivity. In this specific case, the conductivity was first averaged over a range 
6Vg ~ 20 mV that was larger than the average fiuctuations in gate voltage and centred 
around the chosen value for Vg. Minimization procedures were then performed. We did 
not observe any substantial variation of s in this case but a gate voltage independent 
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hopping exponent s ~ 1/3, indicating non- interacting hopping. This clearly shows that 
the crossover observed when varying Vg is progressively disappearing as the ions get 
closer to the Si-Si02 interface. The presence of impurities in the silicon oxide effectively 
creates disorder at the Si-Si02 in the form of potential fluctuations whose amplitude 
is proportional to l/hf and whose length scale varies as hi [19]. Close to the Si-Si02 
interface, these fluctuations are short-range and are pinning strongly the electrons at the 
potential minima. However, this is the long-range character of the Coulomb interaction 
that is responsible of the formation of the soft Coulomb gap in the density of states 
[20] . Thus, hopping is not-interacting and Mott law prevails in this case. The variation 
of s with Vg shows a notable change, especially with the ratio W/U where W is the 
disorder strength and U the on-site energy. The presence of Mott VRH at high disorder 
strength is compatible with W > U characteristic of Anderson localization [21]. But, 
when W < U, a transition to a Mott insulator is expected p2]. So it is legitimate 
to investigate the possibility of an Anderson-like transition in this device [23]. In this 
case, we would expect an increase of the localization length at the transition. Still, the 
determination of from a (T) is difficult because of the variation in s with position of 
the ions in the oxide and the dependence of the density of state at the Fermi level with 
Vg. However, studies at constant temperature and in magnetic field allow getting access 
to the variation of ^ independently. 



3.2. Magnetoconductivity 

In all cases, the magneto-conductivity (MC) is negative up to 10 T, as expected for 
insulating materials [211 EB] and is well described by the general equation 

In a (5) ~ - (B/BoT with = {h/eBof (2) 

The use of the exponent 7 is a convenient way of assessing the non-linearity of the 
MC. Its value ranges from 0.6 to 2 depending on the gate voltage range and position 
of the ions in the oxide. It is interesting to notice that, at specific values for Vg, the 
quadratic behaviour (7 = 2) may extend up to 10 T whereas at others, the linear MC 
(7 = 1) may be present down to OT (Fig. 4). The quadratic term is associated with an 
increase of the localization due to the shrinkage of the electron wave functions centred 
at the impurity site as well as interference effects due to backward electron paths [25]. It 
has already been observed in similar devices where the two-dimensional impurity band 
was formed in the inversion layer [26]. It is classically associated with closed orbits 
and present as long as the hops between the initial and final states are enclosed in a 
cyclotron orbit, e.g. 



e « h'/r (3) 
where /^^ = h/eB and r is the hopping length. 

This condition generally applies at low magnetic field but equivalently in a strong 
localized regime. As for the linear term, it has been observed experimentally by many 
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Figure 4. a (B) at 290mK for Vg=0.7V, 0.58 V and 0.55 V (from top to bottom) 
after DO, Vg=0.6V, 0.47V and 0.33V (from top to bottom) after D2 and Vg=0.46V, 
0.27 V and 0.2 V (from top to bottom) after D4. Best fits are sliown witli dotted lines. 
Dashed lines indicate the noise level. 
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Figure 5. 7 (Vg) and Ib,, (Vg) at 290 mK for DO, D2 and D3. 



authors [27t [28] and explained by the presence of strong random potential fluctuations 
in a two-dimensional electron gas [291 SO] leading to incomplete orbit scattering. In 
the hopping regime, this gives rise to sublinear terms in the MC under appropriate 
conditions ^31j. In our device both phenomena contribute to the MC, and the variation 
of {Vg.W^ reflects the relative strength of the two (Fig. 5). Close to the metal- 
oxide interface, there is a smooth variation from a quasilinear MC for Vg > 0.60 V to a 
quadratic MC for Vg < 0.57 V. At low electron density, Eq. 3 is fulfilled and 7 = 2. At 
higher electron density, the system is no longer strongly localized and ^ increases such 
that it invalidates Eq. 3. The distance between scattering regions diminishes giving rise 
to the linear or sublinear terms in the MC. When the ions are close to the silicon-oxide 
interface, no clear variation of 7 with Vg is observed. Interference effects also contribute 
to the fluctuations in the MC, so that even the use of averaging procedures, as described 
previously, cannot provide a smooth variation for 7 fVgj, but surprisingly, fVgj is 
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linear. At intermediate positions in the oxide, the smooth variation in s, as observed 
for DO, is not present. Instead a minimum is observed in 7 {Vgj in case of drift D2 

for Vg ~ 0.47 V, a value close to the crossover voltage in s i^gj- A similar but weaker 
minimum is also present in the case of the drift D3 at Vg ~ 0.50 V but none for drift 
Dl. In presence of a magnetic field, single hops occur between impurity sites located 
within a cigar-shaped area of dimension r and {r^Y^'^ |32]. Because of the dependence 
of the hopping length on the localization length, the magnetic flux through the hopping 
area is proportional to where a > 0. The magnetic field length as defined in Eq. 
2 is thus essentially dependent on the electron density uq- and the localization length 
[30] so that a non-monotonic variation in Ibq (^) could be attributed to a variation in 

^ {Vgj- In the case of drift DO, Ibq is nearly gate voltage independent indicating that 
the dependence of the hopping length on uq- is weak and so ,^ ~ n^J^ [30j whereas ^ 
is maximum at the crossover in case of drift D2. Such an increase has been associated 
to the presence of an impurity band in low-doped devices [3] but in the present device, 
we did not observe any enhancement of the conductivity. Alternatively, it is possible 
to interpret an increase of ^ by the occurrence of a phase transition. Because of the 
dependencies observed in Figs. 3 and 5, the physical process that is responsible for these 
is likely to rely on both the electron-electron interaction and disorder relative strengths 
[33]. 

4. Conclusion 

We have analysed the electronic properties of sodium-doped silicon MOSFETs by 
transport measurements in both temperature and magnetic field. We have observed 
a complex change in transport mechanism from a Mott hopping to an ES hopping 
regime as a function of the position of the ions in the oxide and the electron density. In 
gate voltage, the variation of the hopping exponent is attributed to the replenishment 
of the Coulomb gap due to the decrease of the localization length with electron density. 
This is possible because of the presence of a long impurity band tails in the density 
of states. This crossover only appears when the ions are within a certain distance to 
the silicon-oxide interface. Close to the metal gate, no crossover is observed and ES 
VRH is present at all sub-threshold gate voltages whereas when the ions are pinned 
at the silicon-oxide interface, the transport is governed by the Mott non-interacting 
hopping. Magnetotransport measurements showed an increase in the localization length 
at the crossover, suggesting the presence of a possible disorder-driven transition in these 
systems, like the Mott-Anderson transition. 
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